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Abstract: We assess the prospect of observing a neutral Higgs boson at hadron 
colliders in its decay to two spin-zero states, a, for a Higgs mass of 90 — 130 GeV, when 
produced in association with a W or Z boson. Such a decay is allowed in extensions 
of the MSSM with CP-violating interactions and in the NMSSM, and can dominate 
Higgs boson final states, thereby evading the LEP constraints on standard Higgs 
boson production. The light spin- zero state decays primarily via a — > bb and t + t~, 
so this signal channel retains features distinct from the main backgrounds. Our study 
shows that at the Tevatron, there may be potential to observe a few events in the 
bbr + T~ or bbbb channels with relatively small background, although this observation 
would be statistically limited. At the LHC, the background problem is more severe, 
but with cross sections and integrated luminosities orders of magnitude larger than 
at the Tevatron, the observation of a Higgs boson in this decay mode would be 
possible. The channel h — > aa — > bbbb would provide a large statistical significance, 
with a signal-to-background ratio on the order of 1 : 2. In these searches, the main 
challenge would be to retain the adequate tagging efficiency of 6's and r's in the low 
Pt region. 
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[51. Summary 



1. Introduction 

The elucidation of the mechanism leading to the origin of mass of all observed elemen- 
tary particles is one of the main goals in high energy physics. The simple Standard 
Model (SM) picture, based on the spontaneous breakdown of the electroweak sym- 
metry by the vacuum expectation value (vev) of an elementary Higgs field, seems to 
lead to a picture that is consistent with all experimental observables, provided the 
Higgs boson mass is smaller than about 250 GeV. Moreover, the best fit to the pre- 
cision electroweak observables measured at the LEP, SLC and Tevatron experiments 
lead to values of the Higgs mass of the order of or smaller than the present bound 
coming from direct searches at LEP, m# SM > 114 GeV [Q. 

In spite of the extraordinary good agreement of the experimental observations 
with the SM predictions, there are many theoretical motivations to go beyond the 
SM description. Several extensions of the SM exist in the literature, and in most of 
them the Higgs sector is extended to a more complicated structure, often including 
at least two Higgs doublets. The requirement of preserving the good agreement with 
experimental data can be easily fulfilled in extensions, like supersymmetry, in which 
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the effect of the additional particles on the precision electroweak observables rapidly 
vanish with increasing values of the new particle masses. An extension of the Higgs 
sector will generally require a revision of the direct and indirect limits on the Higgs 
mass. In particular, the direct search for Higgs bosons may be affected by additional 
decay modes that are beyond the ones analyzed by the LEP collaborations. 

As an example, let us consider the minimal supersymmetric extension of the SM 
(MSSM). In the MSSM, there is an additional Higgs doublet, leading, in the absence 
of CP-violation in the Higgs sector, to two CP-even and one CP-odd Higgs boson 
states. At large values of tan/3, the ratio of of the two Higgs doublets vev's, one 
of the CP-even Higgs bosons acquires SM properties, while the second Higgs boson 
may be produced in association with the CP-odd Higgs boson state. In addition, 
the masses of the non-standard CP-even Higgs and the CP-odd Higgs are close to 
each other. Under these conditions, the mass bound on the SM-like CP-even Higgs 
is similar to the SM one, while the CP-odd and the second CP-even Higgs boson 
mass bound reads rrih > 90 GeV . 

In this paper, we will depart from these simple assumptions, by breaking the 
mass relations that appear in the simplest supersymmetric models, and studying the 
consequences of such modifications of the parameters of the theory. Indeed, while it 
has been a common belief that the Higgs boson will be eventually discovered at the 
upcoming LHC experiments, one would like to fully utilize the potential to search for 
the Higgs bosons at the Tevatron in these non- conventional scenarios as well. Non- 
standard mass relations are already present in extensions of the MSSM including an 
additional singlet (NMSSM [|, |§ and other extensions [§ §, % |, |, [TOf), or when 
explicit CP- violation exist in the Higgs sector flTT| . In these cases, the SM-like Higgs 
(h) may dominantly decay into a pair of lighter Higgs (a), an admixture of CP even 
and odd states with a dominant CP odd component. (The precise fraction of the CP 
even or odd component is not crucial in the present study.) Therefore it is possible 
that the Higgs escaped detection at the LEP experiments by avoiding the usual decay 
modes such as h — > bb, t + t~, WW* and ZZ*, and the lower limit on Higgs mass 
should be re-evaluated. The LEP collaborations have already analyze the possible 
constraints on Higgs boson production arising from this new decay mode |12|, |TjJ. 
We shall use the results of these analyses as a starting point for our study. We are 
interested in analyzing the sensitivity of the Tevatron and the LHC experiments in 
the search for a light, SM-like Higgs boson with such an exotic decay mode. 

We consider the case where the SM-like Higgs boson decays into a pair of spin- 
zero states, h — > aa, which in turn cascade into a heavy fermion pair a — > bb or 
a — > t + t~. These Higgs-to-Higgs decay modes have been studied extensively in 



NMSSM [14, 15, 16] at the LHC, together with even more complicated cascades. 



Most of these studies indeed take advantage of the dominant production modes of 
the Higgs boson at hadron colliders, i.e. gluon fusion and weak boson fusion, but 
encounter large SM backgrounds. We therefore consider the Higgs signal produced 
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in association with a W or Z boson, where the leptonic decays of the weak bosons 
will provide a clean trigger, and will significantly reduce the background as well. 1 

The significance of associated production was also stressed in Ref. for h — > aa 
decays in NMSSM at the LHC For similar parameter choices, our results lead to a 
comparable or slightly better reach for Higgs boson searches to the ones obtained 
in earlier studies |TJ|, [T5], [Tj| on gluon fusion and weak boson fusion productions at 
the LHC, while showing more discovery potential at the Tevatron than in previous 
studies. 



2. Signal Processes and Parameter Choices 
2.1 Signal Processes 

For the two spin-zero states, a's, the combinations of decay products we can search for 
are Ab, 262r and At. The At mode is usually suppressed by the branching fractions, 
unless m a is below the bb decay threshold p0| . We thus concentrate on the two 



channels 262r and Ab next. The signal events being searched are 

Wh ^ lvh0a ^ _ (2.1) 

L lv u bb,T + T 

7 , ,+ ,- (l + l~,bb,bb 

Zh^m >a a->| i+I _ SjT+T _ (2.2) 

with I — e,fi. The channel Z —>■ vv decays into neutrino pairs can also be considered, 
while the triggering could be large missing energy, plus r's or 6's. 

2.2 Parameter Choices 

We would like to perform a relatively model-independent search for the signal, there- 
fore the Higgs masses, branching fractions and couplings to the weak bosons are 
employed as input parameters. Direct searches for a Higgs boson with SM-like cou- 
plings to the gauge bosons, in a model and decay mode-independent way, leads to 
a lower bound on m h of about 82 GeV |21| with full SM coupling to Z. On the 



other hand, the proposed search is expected to become inefficient for m h > 130 GeV, 
since the standard decays into the WW* and ZZ* channels are still expected to be 
dominant. Therefore, the optimal setting to detect the Higgs decaying into an aa 
pair is to have the mass within the range of 90 — 130 GeV. The choice for m a can 
be more flexible. As long as rrih > 2m a and m a > 2m& to kinematically allow the 
decays h —>■ aa and a — > bb, our methods are rather insensitive to the mass choices. 



1 An initial analysis by us at the Tevatron was reported earlier in Ref. |Q. While this current 
work was in process, another similar analysis for the 46 channel at the LHC appeared [jl8| . For the 
overlap with that work at the LHC, our studies included more background analyses, more realistic 
6-tagging effects, and a broader parameter scan. 
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In a generic model, the Wh/Zh production rate differs from that in the SM. The 
change can be characterized by a prefactor nf lWW (k\ zz ), where Khvv is the coupling 
strength of Higgs to vector boson V relative to that in the SM. The production cross 
section can thus be written in terms of the SM result with an overall factor to account 
for the modification of the coupling 

a(Vh) = n 2 hVV a SM (Vh). (2.3) 

We are interested in the range of k 2 ~ 0.5 — 1.0, so that this Higgs contributes to 
the electroweak symmetry breaking and consequently the associated productions are 
still sizable. 

In order for the h — > aa decay to be dominant and thus escape the LEP bound, 
BR(h — > aa) is required to be close to unity. For instance, in the NMSSM, BR(h — > 
aa) > 0.9 turns out to be very general in terms of the naturalness of c in the trilinear 
coupling term (cv/2)haa [p2j1 . Moreover, if the down quark and lepton coupling to 
the Higgs is proportional to their masses, then BR(a — > bb) and BR(a — * t + t~) 
are set to be 0.92 and 0.08, respectively. In general, however, the relations between 
the coupling and the masses may be modified by radiative corrections, which can 
lead to a large increase of the BR{h — > rr) p3[ . The representative values and the 
ranges of the parameters are summarized in Table [1], all allowed by constraints from 
LEP ||i"2"| , T3\, except for the region near rrih ~ 90 GeV when both a's are assumed to 
decay into two bottom quarks. Parametric consistency with the LEP results is also 
discussed in detail within the NMSSM framework. 







representative 


considered 




parameters 


value 


range 


masses 


m h 
m a 


120 
30 


90-130 
20-60 


coupling 


K hVV 


0.7 


0.5-1.0 


branching 
fractions 


BR(h — > aa) 
BR(a bb) 


0.85 
0.92 


0.8-1.0 
0.95-0.50 


BR(a -> t + t-) 


0.08 


0.05-0.50 


262r channel 




0.088 


0.038-0.50 


46 channel 




0.50 


0.10-0.90 



Table 1: Parameter choices for h 
section. 



aa decays. The C 2 factor is defined in the next 



3. h — > aa at the Tevatron 
3.1 The 262r Channel 

Including the decay branching fractions for aa — > bbr + T~, we obtain the cross section 
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as 

^2b2r = cr(Vh) BR(V) 2BR(h -> aa)BR(a -> bb)BR(a -> r + r~). (3.1) 

where = 0.213 (0.067) is the leptonic branching fraction of (Z) decay into 

/ = e, /i. 

The overall factor modifying the SM result in Eq. ( |3.1|) , 



C 2 2 62t = 2K 2 hvv BR(h -> aa)BR(a -> bb)BR(a -> r + T~), (3.2) 

corresponds to the pro cess- dependent C 2 factor defined in the DELPHI search [|T2j|, 
and the 6*95 factor in the comprehensive LEP analysis [jnj. Our parameter choice 
(range), as listed in Table [TJ, is equivalent to a C| 62r of 0.088 (0.038—0.50), consistent 
with the bounds for a large range of our mh, m a choices set forth in Refs. |12| , |13 |. A 



value of 0.088 for C\ h2r is assumed for all numerical evaluations from here on, unless 
explicitly noted otherwise. 

3.2 Signal Event Rate for the 262r Channel 

The associated production of pp — > Wh usually features a larger cross section than 
that of Zh, and the leptonic branching fraction of W is about 3 times larger than 
Z's. For illustration purposes, we choose to present our detailed studies for the Wh 
channel henceforth, although we will include the Zh channel in our results. 

Total cross sections for pp — > Wh and pp — > Zh at hadron colliders have been 



calculated with QCD and electroweak corrections included [^, [26|, |27| in the SM. 
Hence we get 

a SM (Wh) BR(W -> lv{) ~ 85 (24) fb (3.3) 

at y/s = 1.96 TeV for m h = 90 (130) GeV. 

Including the branching fractions and couplings, the cross section of the signal 
in Eq. (Q is 

cr 2 b2r ~ 7.5 (2.1) fb for C 2 = 0.088 (3.4) 

as illustrated in Fig. [1]. The solid curve on top represents the total cross section 
for Vh production, with V decaying leptonically, but without any cuts. The dashed 
curve represents the cross section after adjusting for the couplings and branching 



fractions, as in Eq. (|3.1| ). Cross sections for Zh are also plotted for completeness. 
3.3 Background and Cuts for the 2b2r Channel 

The main advantage for considering the process Wh is the possible background 
suppression due to the clean final state from the W leptonic decay: an isolated 
charged lepton (I = e, /1) plus large missing transverse energy. We thus require the 
following initial acceptance cuts at the Tevatron [ f28]j 

p T (/)>15GeV, |??(/)|<2.0, $ T > 15 GeV. (3.5) 
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Wh -> lvbbx X at Tevatron 
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Figure 1: Cross sections of Higgs signal at the Tevatron in the 2b2r channel produced 
by Higgs-strahlung with a leptonically decaying W (left) or Z (right). The four lines 
from the top to the bottom correspond to, respectively, the SM cross section (solid line); 
adjusted for C\ h2r = 0.088 and m a = 30 GeV (dotted line); including the acceptance cuts 
Eqs. Q3.5|) — ( |3.8D (line with circles); and further including tagging efficiencies Eq. Q3.9| ) (line 
with crosses). The shaded bands correspond to variations of the final results (line with 
crosses) for values of C 2 within the range considered in Table |l] and taking into account 



the LEP constraints [13] 



The events have yet to further pass the acceptance cuts, or to have the taus and 6's 
tagged. Both help suppress the SM backgrounds, while bringing significant reduc- 
tions to the event rate as well. Our challenges are to retain as many signal events 
as possible, and to control the backgrounds from various sources. Throughout the 



paper, we adopt the Monte Carlo program MadEvent j29[ for our background simu- 
lations at the parton level. 



b and r Tagging We wish to identify events with 5 particles plus missing energy 
in the final states: bbr + r~lvi. With neutrinos in the decay products, tau momenta 
cannot be fully reconstructed. Therefore we cannot reconstruct the invariant masses 
m(2r) or ~ m(262r). Instead, the signal should appear as a peak in the m(2b) 
plot, around the value of m a . 

For the jets and other soft leptons in the events, the following basic cuts are 
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employed to mimic the CDF |30| detector acceptance, for jets [[Hp 



PrO-)>10GeV, < 3.0, (3.6) 

and for r-candidates |3^| : 



p T > 10, 8,5 GeV for r^,r e ,r M , \r)\ < 1.5. (3.7) 



7") 



where r e , r M and stand for the visible decay products of r — > ev e v T , r — ► [w^v, 
and r — > hadrons + i/ r , respectively, and an isolation cut 

Ai? > 0.4 (3.8) 

between leptons, r's and 6-jets. After the acceptance cuts, 10 — 25% of the signal 
events survive, and the cross section becomes 0.85 (0.57) fb for = 90 (130) GeV 
with the given set of input parameters [C 2 ~ 0.088). The cross sections passing 
acceptance are plotted in Fig. [I] versus the Higgs mass, represented by the circled 
curve. At this level, the cross section is below 1 fb. 

The 6- and hadronic r-tagging efficiencies and the kinematics are taken to be 

e b = 50% for E^ 1 > 15 GeV and \ij jet \ < 1.0 , 

e T = 40% for E vis > 20 GeV and \r}\ < 1.5 . (3.9) 



Outside these kinematical regions, the tagging efficiencies drop off sharply |33|, 
We decide to tag one b and one tau. The energies for a jet and a lepton are smeared 
according a Gaussian distribution. The energy resolutions are taken to be 

^- re « ffi 5%. ^ = ±^®1%. (3.10) 



Ej sJE, El iffi 

The missing energy is reconstructed according to the smeared observed particles. No 



further detector effects are included |35]. 



Irreducible Background The dominant source of the irreducible background, 
with the same final state as the signal, 

W Z*/t*(-> t + t-) 66, (3.11) 

has the bb pair from a virtual gluon splitting, the t + t~ pair from an intermediate 
Z*/'-/* and the charged lepton plus missing energy from a W boson. Our simulations 
show that the largest contribution come from events with the Z* almost on-shell, 
while the r + r~ pair from a virtual photon can be more easily confused with the 
signal. After applying the acceptance cuts, the irreducible background is estimated 
to be around 0.01 fb, which is very small compared to the signal size. It is essentially 
absent given the luminosity expected at the Tevatron. 
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Reducible Background Reducible background arise from jets mis-identified as 
b's, or as hadronically decaying taus. The mistag rate from a light quark is taken to 
be 0.5 — 1.0% for tau and 0.5% for b, respectively |33|, A charm quark has higher 



mistag probability to fake a b quark, that we take to be 10% |36[. In addition, the 
experiments cannot distinguish directly produced electrons (muons) from leptonically 
decaying taus. Thus the reducible backgrounds considered in our study are listed 
below. 

• The background due to misidentified bottom comes from the process 2r2j I + 
$ T , which has a cross section of 5 fb. Considering the mistag rate and the 
additional cuts, it contributes 0.02 fb to the background events. 

• The background due to misidentified r differs for different decay modes of r's: 

— For TiThlbl$ T (2l2brh$T)) it comes from 2r2bj with $ T from the leptonic 
decays of both taus. The contribution is estimated at 0.003 fb. 

— For ThTh2bl$ T) the background comes from 2j2bl$ T and is estimated at 
30 fb. It is then reduced by the tau-mistag rate, the b-tagging rate, 
and their associated cuts. In the continuum distribution of m(2b), it 
is below the level of the resonant signal. Within the mass window of 
lOGeV < m{2b) < 70 GeV, this background accumulates to 0.04 - 0.09 
fb, depending on the r mistag rate considered. 

• The backgrounds from both a mistagged tau and a mistagged b mostly come 
from the 4j7$ T events, which has a cross section of about 16 pb. After the cuts 
and folding in the mistag rates, this contributes 0.03 — 0.05 fb of background 
in 10 GeV < m{2b) < 70 GeV, depending on the r mistag rate considered. 

The two bottom-quarks in the final state coming from the Higgs boson decays 
should have an invariant mass equal m a . If enough data were available, one would 
be able to observe an excess of events in the m(2b) mass distribution. However, 
this procedure is heavily limited by statistics. For instance, with a window cut of 
m a ± 10 GeV on m(2b), the reducible background can be a factor of 3 to 5 smaller 
than the signal, but unfortunately, the cuts and the tagging efficiencies together 
reduce the signal greatly to about 0.11 fb for Wh and 0.05 — 0.07 fb for Zh, with 
C 2 ~ 0.088 as shown in Fig. [I] by the crossed curve. The shaded band represents the 
range of parameters allowed by our choice of C 2 ~ 0.038 — 0.50, consistent with the 
LEP constraints. With an optimistic value of C 2 ~ 0.50, the cross section is 0.68 
fb, and we would expect to see about a couple of signal events with an integrated 
luminosity of a few fb -1 . 

To illustrate a most optimistic situation in terms of kinematical considerations, 
we explore the optimization between m a and rrih to obtain the largest signal rate. 



- 8 - 



Wh -> lvbbbb at Tevatron 
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— Standard Model 

— Adjusted for K.BRs 
G-O Acceptance 

>^< Tagging 




110 120 
ra (GeV) 



100 110 

m h (GeV) 



130 



Figure 2: Cross sections of Higgs signal at the Tevatron in the 4b channel produced by 
Higgs-strahlung with a leptonically decaying W (left) or Z (right). The four curves from 
top to bottom correspond to, respectively, the SM cross section (solid line); adjusted for 
C 2 h = 0.50 and m a = 30 GeV (dotted line); including the acceptance cuts Eqs. (j3li|), fl3.6| ) 
and ( |3.8| ) (line with circles); and further including tagging efficiencies Eq. ( |3.9| ) (line with 
crosses). The shaded bands correspond to variations of the final results (line with crosses) 
for values of C 2 within the range considered in Table [l] and taking into account the LEP 



constraints [13| 



The signal loss is mainly due to the softness of the b and r's, therefore most events 
are rejected from the lower p? threshold. Increasing m a would stretch the pt dis- 
tributions to the higher pt end. To achieve this without significantly affecting the 
decay phase space of h, we set 



(m h - 10GeV)/2, 



(3.12) 



which resulted in more than doubling the signal rate with respect to the m a = 30 GeV 
case as our default presentation throughout. In this case the signal cross section is 
~ 0.28 fb for C 2 = 0.088, and ~ 1.6 fb for C 2 = 0.50, which is still challenging for 
observation with the Tevatron's projected luminosity. 

3.4 The 46 Channel 

As we mentioned earlier, the light spin-zero state a mostly decays into bb (50 — 95%) 
or r + r~ (5 — 50%). The t + t~ channel can be dominant when a is very light, i.e. 
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m a < 2m&. However it would be difficult to observe Higgs in the 4r mode, first 
because of the increasing background near the lower end of m(rr), and because of 
the difficulty in resolving highly-collimated tau pairs. These scenarios involving a 
very light a are among the difficult ones for NMSSM Higgs discovery discussed in 
Ref. \TE\. A relevant 4r study at the Tevatron under such scenario can be found 



in Ref. J20]. A very light a could also be probed through Upsilon or even J/\l/ 
decays [57]. For the m a > 2m b case, we will next look for the Higgs in the 4b 
channel. 

Similar to the 262r mode, the 4b cross section is 

a Ab = a (Vh) BR(V) BR(h aa)BR(a -> bb) 2 , (3.13) 
from which we extract the C 2 factor 

C\ h = K 2 hVV BR(h -> aa)BR(a -> bb) 2 . (3.14) 

The 4b mode is usually enhanced by the large branching fractions of the decay 
of a into bottom quarks. The ratio C\ h jC\ mr ranges in 9.5 — 0.5 for BR(a — > 
tt) ~ 0.05 — 0.50. The value of Cf b itself does not vary greatly with the branching 
fractions that are obtained within our choice of parameters, Table [l]. Despite larger 
background for this mode than for the 262r mode, the enhanced rate suggests this 
to be a more viable mode. 

The parameter choices for the 4b channel are also given in Table. [1]. Running 
parallel to the 262r channel, we plot the cross sections in Fig. |2|. The shaded bands 
show the LEP constraints disfavoring the lower end of the rrih range. We find the 
signal rate after acceptance cuts to be 10.7 — 4.7 fb (the circled curve) for rrih = 90 — 
130 GeV with C 2 ~ 0.5. After tagging three bottom jets (for reasons explained below) 
and imposing appropriate additional cuts, the cross section becomes 0.54 — 0.38 fb 
(the crossed curve) for rrih = 90 — 130 GeV. 

We again adopt the basic acceptance cuts and the 6-tagging requirements as in 
the previous section. The background for this mode arises from 4bl$ T , 3bjl$ T ,2b2jl$ T , 
b?>jl$ T and 4jl$ T events. For the four b's in our signal, tagging two will not be suffi- 
cient, as background from 2b2jl$ T events can fake the signal without any mistagging 
involved. Therefore we demand that at least three bottom jets be tagged. 

The irreducible background 4bl$ T) though much larger than that in the 262r 
mode, is still manageable. We find the cross section to be 0.23 fb after basic ac- 
ceptance cuts. Like the signal events, it suffers similar reductions from tagging and 
further cuts, which brings it down to 0.02 fb. With tagging for 36's, the 3bjl$ T 
events cannot be effectively distinguished from the signal either. They contribute 
about 0.003 fb to the background. 

The reducible backgrounds from 2b2jlp T and 4jl$ T events have the same sources 
as that in the 262r mode and the mistag rates of a light jet to b and to r are 
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comparable. The tagging on the 3 rd 6 brings this background down significantly. In 
total, they contribute about 0.07 fb to the background. Another background source 
is from 2b2cl$ Tl which is approximately 2.5 times as large as that of 4bl$ T at the 
Tevatron. With a 10% mistage rate and after the acceptance cuts, it is reduced to 
0.007 fb. Finally 3bjl$ T and b3jl$ T backgrounds combine to contribute less than 
0.003 fb. 

Having tagged three of the four bottom quarks, we identify the fourth bottom 
as the hardest untagged jet in the event. We expect the signal to appear as a 
peak in the invariant mass m(b\, 62) and 771(63, 64) distribution. However, pairing the 
four b jets can be complicated due to combinatorics. We assign the two pairs by 
minimizing their mass difference 777,(61,62) ~ 777(63,64) and record both values each 
with a half weight. We present the signal versus the background distributions of 
the reconstructed masses and m a in Fig. [| as the invariant masses of four 6-jets 
and of two 6-jets. With a simple cut on the 777(46) invariant mass, 777(46) < 160 GeV, 
dictated by our search for a light Higgs boson with mass smaller than about 130 GeV, 
the overall signal to background ratio can be about 10 with C 2 = 0.50, 777 a = 30 GeV 
and 777^ = 90 - 130 GeV. 

To summarize our study at the Tevatron, we claim that the signal channels of 
Eqs. ( [2.1D and ( |2.2| ) have distinctive kinematical features (see Fig. ||) with negligible 
SM backgrounds and the signal observation is total statistically dominated. For the 
262r mode, one can reach a cross section of about 0.05—0.7 fb as shown in Fig. [I], 
while for the 46 mode, we have the cross section in the range of 0.1—1 fb as shown 
in Fig. |2|. If the h and a masses happen to be related in an optimal way (Eq. (|3 . 1 2|) ) 
we can gain an increase in the signal rate by a factor of ~ 1.8 and 2.5 for the 46 and 
262r channels. 



4. h — > aa at the LHC 

At the LHC, weak boson-associated Higgs production rate is about 10 — 15 times 
that at the Tevatron in the mass region we are interested in. With the same C 2 
factor, signal events passing through acceptance also take on this ratio. The (QCD) 
background, on the other hand, can be 100 times larger than at the Tevatron. This 
requires a substantial jet rejection rate. These cross sections are plotted in Fig. f|. 
Cuts on the triggering leptons and/or missing energy are taken to be 

Pr(Z)>20GeV, rj(l) | < 2.5, $ T >20GeV. (4.1) 



The following cuts and efficiencies for tagging are assumed [38 



e b = 50% for E J T et > 15 GeV and \rj jet \ < 2.0 , 

e r = 40% for E vis > 15 GeV and < 2.5 . (4.2) 
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Figure 3: Higgs signal (double-hatched) on top of the sum of the backgrounds at the 
Tevatron in the 46 decay channel together with a leptonically decaying W. The invariant 
mass of four (left) and two (right) 6-jets are shown. Values of C\ h = 0.50, = 120 GeV and 
m a = 30 GeV are understood. From bottom to top, the background histograms indicate 
the accumulative sum of 2b2cW + 262 jW + 361 jW and 2b2cW + 262 jW + SbljW + 4bW. 

The jet rejection rate is better than 1/150 for tagging a b or a r, except in the 15 — 
30 GeV pt range where it is taken to be ~ 1/30, as there exists strong tension between 
tagging efficiencies and the jet rejection rates, especially near the low p? range. Note 
that the jet rejection rate will only be accurately known after understanding the 
detectors with examining the real data. Again, in our simulations, the energies for a 
jet and a lepton are smeared with the Gaussian resolutions 

AE,j 50% M AE t 10% M , x 

J ~ © 3%, — - = —. = © 0.7% . (4.3) 



Ej JE 3 E x y[E x 

The missing energy is reconstructed accordingly. 
4.1 The 262r Channel 

Similar to the Tevatron case, the irreducible background of Eq. ( |3.11| ) is small after 
the acceptance cuts and the tagging requirements, contributing only 0.07 fb. The 
reducible background, however, poses a much more severe problem at the LHC For 
example, the 2b2jl$ T events are estimated to be around 11 pb, compared to 50 fb at 
the Tevatron. Thus for the 262r mode, a jet rejection rate of 1/150 would give rise 
to a background of 92 fb, compared to the signal size of about 1 fb (or up to ~7 fb 
when maximizing C 2 ). The 4jl$ T events also contribute 43 fb to the background in 
this channel. 
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Figure 4: Cross sections of Higgs signal at the LHC in the 262r (left) and 46 (right) 
channels produced by Higgs-strahlung with a leptonically decaying W. The four curves 
from top to bottom correspond to, respectively, the SM cross section (solid line); adjusted 
for C| 62r = 0.088, C 2 b = 0.50 and m a = 30 GeV (dotted line); including the acceptance 
cuts Eq. ( |4.1| ) (line with circles); and further including tagging efficiencies Eq. Q4.2j ) (line 
with crosses). The shaded bands correspond to variations of the final results (line with 
crosses) for values of C 2 within the range considered in Table |l] and taking into account 



the LEP constraints [13] 



We carry out the analysis similar to the Tevatron case and arrive at a S/B ratio 
of 0.03, with a total signal size of less than 1 fb for = 120 GeV, m a = 30 GeV 
and C 2 = 0.088. The small S/B ratio would require precise control of the systematic 
errors. It can be further improved by tagging one more b or r, at the expense of 
losing up to half of the signal rate. Due to the difficulty of finding a signal in this 
channel, we are led to consider the more promising channel of 46's, where, as we did 
in the Tevatron case, we employ an additional tagging, while still retaining a higher 
signal rate. 



4.2 The 46 Channel 

With a much higher luminosity than the Tevatron and larger cross sections, LHC 
could produce 60 (10 fb -1 ) to over a thousand (300 fb -1 ) Higgs events in the 4bl$ T 
decay channel, assuming a typical C 2 value (C| b = 0.50), as shown in Fig. |. The 
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Figure 5: Higgs signal (double-hatched) on top of the sum of the backgrounds at the LHC 
in the 46 decay channel together with a leptonically decaying W. The invariant mass of 
four (left) and two (right) 6-jets are shown. Constraints of 60GeV < m(4b) < 160 GeV and 
lOGeV < m(26) < 70 GeV are implemented in both plots. C\ h = 0.50, m h = 120 GeV and 
m a = 30 GeV are understood. From bottom to top, the background histograms indicate 
the accumulative sum of 2b2cW, 2b2cW + 2b2jW, 2b2cW + 2b2jW + SbljW, and 2b2cW + 
2b2jW + SbljW + 4bW, respectively. 

46 channel is thus more optimistic for observing the Higgs, even though the back- 
ground still dominates the signal, and the irreducible 4bl$ T background becomes 
non- negligible. We require tagging three of the 6 jets, which would essentially elim- 
inates backgrounds from 4jlp T , and reduces the 2b2jl$ T and lb3jlp T background 
significantly. With three tagged 6-jets, the signal rate is about 5.7 fb (or up to 10 fb 
when maximizing C 2 ). The irreducible background 4bl$ T is 25 fb. The 3bjl$ T 
background is about 16 fb. The reducible background from 2b2jl$ T events is about 
80 fb, and 2b2cl$ T is about 4 fb with a 10% mistag rate for c — > b [[38]. The 4j/j5 T 
background is no larger than 0.2 fb. 

We again present the reconstructed mass distribution for the signal and back- 
grounds in two plots in Fig. [5[ The left and right plots show the invariant mass 
distributions of the 46 and 26 system, where the signal peaks near rrih = 120 GeV 
and m a = 30 GeV, respectively, each with a width less than 10 GeV due to detector 
energy resolution. Similar to the Tevatron case, we assign the two 66 pairs by min- 
imizing their mass difference 771(61,62) ~ 771(63,64) and plot these two masses, each 
with a half weight. 

The dominant 2b2jlEj- background comes from ti production. For ti events, the 
262j system contains all the decay products of a top-quark. Therefore, these events 
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may be efficiently rejected with an upper cut on the m(4b) invariant mass lower than 
the top quark mass, m(46) ^160 GeV, which will not affect the signal we consider if 
the Higgs boson mass is in the region m < 130 GeV. Given our considered range of 
choices, we implement the following constraints in the two distributions: 

10 GeV < m(26) < 70 GeV, 
60 GeV < m(46) < 160 GeV. 

While the former affects the m(46) distribution minimally, the latter reduces the 
background in m(2b) distribution by about 40%. 

Overall, selecting events with these invariant mass constraints, the value of S/ B 
is roughly 1/5 for C\ h = 0.50. Assuming a good understanding of the background, 
one can get an estimate of the statistical significance of the signal. For the rate 
quoted above we obtain a significance, S/yB, of over 3.5a for 10 fb _1 and over 5a 
for 30 fb _1 , as indicated in Fig. [5| If one selects events only in the expected signal 
region, we obtain a S/B ~ 0.41 in the range 100 GeV < m(46) < 140 GeV from 
the m(46) distribution, and a S/B ~ 0.40 in the range 20 GeV < m(2b) < 40 GeV 
from the m(2b) distribution, equivalent to a reduction by about a factor of two of 
the luminosity necessary to achieve the same statistical significances. The challenge 
is for us to understand the background well enough, and to control the systematic 
errors. 

It may be a challenge at the LHC to retain the high 6-tagging efficiency at 
Pt ~ 15 GeV adopted in the current analysis. If a 30 GeV cut on the tagged jets is 
implemented instead, the signal is reduced to 22%, while the background drops to 
about 37% of the values given above. In such case a 3a (5a) signal would require 
an integrated luminosity of around 30 (80) fb -1 . Therefore a good understanding of 
6-tagging efficiencies at low px will be necessary to be able to discover a Higgs in the 
46 channels in the first years of the LHC. 

Before closing this section, a remark is in order for comparing our results with 
a recent similar analysis for the 46 channel at the LHC ||18|| . Their conclusions 
are somewhat more optimistic, largely due to a significantly higher 6-tagging effi- 
ciency assumed (70%). They did not consider the QCD backgrounds of W2b2c and 
W2b2j, which are sub-leading. On the other hand, we neglected the background tibb 



considered in [18| since with the additional energetic W from the top-quark decay 



Ew ~ ^r\A — Mw/rriti this background can be efficiently removed by vetoing the 
extra jet or charged lepton activities from the W decay. 

5. Summary 

The search for a Higgs boson with couplings to the gauge bosons of the order of the 
SM one, and decaying into two lighter CP-odd Higgs bosons states may be performed 
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at hadron colliders for the associate production of Wh, Zh with h —>■ aa (262r or 46). 
The cross sections scale proportionally to C 2 , a factor determined by the product of 
the relevant branching fractions times the ratio of the Higgs production cross section 
to the SM one. Maximal event rates of the two channels are given by different 
values of BR(a — > rr). SM-like aff couplings tend to give small BR(a — > rr), 
thus suppressing the 262r channel and enhancing the 46 channel. In models where 
BR(a — > rr) is large, the 262r channel yields an event rate comparable to the 46 
channel. 

We analyzed the Wh channel in the mass range 90 < rrih < 130 GeV in detail. 
We found that at the Tevatron 

• With only basic cuts, the signal size is 0.7 fb for the 262r channel for C 2 b2r ~ 
0.088 with a negligible irreducible background, and 5 — 10 fb for the 46 channel 
for C\ h ~ 0.50 with a comparable background. With favorable couplings and 
branching fractions, the C 2 factor can be as large as 0.50 for the 262r mode, 
and 0.90 for the 46 mode, and the signal rate is enhanced proportionally. 

• Further cuts and the tagging of 6 and r, necessary to remove the much larger 
reducible background, worsen the signal event rate to around 0.11 fb for the 
262r mode and 0.5 fb for the 46 mode, as summarized in Figs. [I] and §. However, 
the kinematics of the mass reconstruction of m a and rrih can be very distinctive, 
as seen in Fig. |3] for the 46 mode with small background and a couple of total 
events. 

• We also consider the most favorable relations between rrih an d the CP-odd 
Higgs mass m a , which can enhance the signal rate by a factor of 2.5 for the 
262r mode leading to a cross section as large as 1.6 fb with C\ h2r = 0.5, and 
by a factor of 1.8 for the 46 mode leading to a value 1.8 fb for C 2 b = 0.9. 

• There can be another improvement of 15 — 30% by combining Wh events with 
the Zh events, where both Z —>■ 11 and Z — > vv can be included, leading to a 
possible observation of a few events in either 262r or 46 channel, for a Tevatron 
luminosity of the order of a few fb -1 . 

Overall, the signal observation becomes statistically limited. Our study has been 
based on parameters of the CDF detector. One expects the signal observability to 
be enhanced accordingly if results from the DO detector were combined. 

At the LHC, the signal rate increases by a factor of 10, and the background 
increases by two orders of magnitude, compared to the Tevatron. We found that 

• Statistics limitation is no longer a major issue. In the 46 channel alone the 
signal rate is 5.7 fb , and we can easily obtain a signal significance S/y/B 
greater than 3.5 with an integrated luminosity of 10 fb -1 , and over 10 with 
100 fb" 1 . 
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• Similar to the Tevatron study, with favorable couplings and branching fractions, 
the signal rate can be enhanced to be as large as 10 fb with C 2 = 0.9, as seen 
in Fig. |], and S/B can be improved accordingly. 

• The kinematics of the mass reconstruction of m a and trih can be very distinctive, 
as seen in Fig. |5] for the 46 mode, yielding a statistically significant signal. 

• Favorable rrih and m a relations, combinations of the Wh and Zh signals, com- 
binations of the 262r channel with the 4b channel, could all improve the signal 
rate and enhance the potential to the eventual discovery of the Higgs boson. 

The main challenge would be to retain the adequate tagging efficiency of b's and r's 
in the low region. 

We point out that our background analysis is based on the leading order partonic 
calculations in MadEvent. More accurate estimate of the background distributions 
would be important to claim a signal observation. More realistic simulations includ- 
ing the detector effects are needed to draw more convincing conclusions. 
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